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' Abstract 

^ , We investigate an additive perturbation of a complex Wishart random ma- 

trix and prove that a large deviation principle holds for the spectral measures. 
The rate function is associated to a vector equilibrium problem coming from 
' logarithmic potential theory, which in our case is a quadratic map involving 

the logarithmic energies, or Voiculescu's entropies, of two measures in the pres- 
ence of an external field and an upper constraint. The proof is based on a 
, two type particles Coulomb gas representation for the eigenvalue distribution, 

p ^ ' which gives a new insight on why such variational problems should describe the 

limiting spectral distribution. This representation is available because of a Nik- 
ishin structure satisfied by the weights of the multiple orthogonal polynomials 



hidden in the background. 



1 Introduction and statement of the results 
1.1 Introduction 



, The study of the large deviations for the spectral measures of large random matrices 

has started with the work [6] of Ben Arous and Guionnet, and continued with many 
■ extensions, see e.g. [7, 34, 26, 12, 31], which now cover all the so-called unitary 

! invariant matrix models, and actually the larger class of /3-ensembles. The proof 

of such large deviation principles (LDPs) is based on the fact that an explicit and 
tractable expression is available for the joint eigenvalue distributions, which is a 
consequence of the unitary invariance. A common feature shared by these random 
matrix ensembles is that the rate functions governing such LDPs, which are maps 



on the space of probability measures, are given by the logarithmic energy functional 

log] — - — 7dn{x)dfi{y), (1.1) 

F - y\ 

plus a linear term in the probability measure. The latter functional (1.1) is the main 
object of study in logarithmic potential theory, and has moreover been interpreted 
up to a sign by Voiculescu as the free entropy, a free probability equivalent of the 
Shannon's entropy in classical probability [43], see also [9, 34, 33]. 
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More recently, much attention has been given to perturbed matrix models where 
one has broken the unitary invariance by the addition, or multipHcation, of an ex- 
ternal deterministic matrix, and also multi-matrix models. It is a highly non-trivial 
problem to establish in full generality that a LDP still holds for such matrix mod- 
els, because of the complex dependence between the eigenvalues and eigenvectors. 
By developing an appropriate non-commutative Ito calculus, Cabanal-Duvillard and 
Guionnet obtained a LDP upper bound for the spectral measures of a large class of 
matrix valued stochastic processes [14]. It has been later extended to a full LDP 
by Guionnet and Zeitouni [29, 30], and a LDP for perturbed or multi-matrix mod- 
els actually follows by contraction principle. The price to pay for such a level of 
generality is a quite complicated rate function, but it is worth mentioning that it 
is known to reduce to the logarithmic energy in the unitary invariant case (i.e null 
perturbation), see [15, Section 5.1]. 

In this work, we shall follow a different path and explore the large deviations of a 
perturbed matrix model through its connection to multiple orthogonal polynomials 
(MOPs). Indeed, while the unitary invariant matrix models are known to be related 
to orthogonal polynomials [35], it has been observed by Bleher and Kuijlaars that 
perturbed matrix models benefit from a connection with MOPs [11], in the sense 
that the average characteristic polynomial of the random matrix is a MOP with 
respect to appropriate weights and multi-index. Such relation also holds for multi- 
matrix models [24], see also [37] for a survey. On the other hand, the limiting zero 
distribution of certain classes of MOPs can be described in terms of the solution of a 
vector equilibrium problem [3, 40] : given d > 1 and a dx d real symmetric positive 
definite matrix C = [ctj], minimize the functional given by 



plus linear terms in (/xi, . . . , fi^), when the vector of measures (^ui, . . . , fid) runs over 
A^mi(Ai) X ••• X Aimai^d), or in some subset thereof. Here A^m(A) stands for 
the set of Borel measures on A C C with total mass m. For a general treatment 
concerning vector equilibrium problems, see [5, 32]. 

A natural question is then to seek if the functionals associated to vector equilib- 
rium problems should be involved as large deviations rate functions. It is the aim 
of this work to answer affirmatively for a particular example that we present now. 

1.2 Non-centered Wishart random matrix 

The model we investigate here is a non-centered Wishart random matrix, which is 
an additive perturbation of the usual Wishart model. Namely, let X = [Xij] be a 
M X N complex matrix filled with i.i.d (non-centered) complex Gaussian random 
entries Xij ~ Mci^ij , 1 / y/N) , where A = [Aij] is a given deterministic M x N 
complex matrix. One can equivalently endow the space A^M,iv(C) oi MxN complex 
matrices with the probability distribution 





(1.2) 
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where Zm,n is a normalization constant and dX stands for the Lebesgue measure 
on A4m,n{C) — M^^"^^. Without loss of generality, A can be chosen in its singular 
value decomposition form. Note that, if Un{C) stands for the unitary group of C^, 
Pat is not invariant under the transformations X i— t- UXV* for given U G UMiC), 
V e Un{C), except if A = 0. 

We are interested in the convergence and deviations of the spectral measure 



N 



1 ^ 



N 
1=1 

where the the eigenvalues of the non-centered Wishart matrix X*X (or 

equivalently the squared singular values of X) with X drawn according to Pjy. It 
is a random variable taking its values in A^i(R+), that we equip with its weak 
topology. 

This matrix model has been extensively studied in the statistic and signal pro- 
cessing literature (see e.g. [41] and references therein), and Dozier and Silverstein 
described the limiting eigenvalue distribution for a large class of perturbations A by 
means of a fixed point equation for its Cauchy-Stieltjes transform [21, 22]. Alter- 
natively, the limiting eigenvalue distribution can be characterized in terms of the 
rectangular free convolution introduced by Benaych- Georges [8]. On the other hand, 
the non-centered Wishart matrix model does not belong to the class of random ma- 
trices for which Guionnet and Zeitouni were able to extend the LDP upper bound of 
Cabanal-Duvillard and Guionnet into a full LDP for the spectral measures {fi^)]\f, 
and to prove such a LDP is in fact still an open problem. 

In this work, we shall restrict our investigation to a particular case and assume 
that M = N + a where a is a non-negative integer, and consider for a > the 
particular type of perturbation 



A 



G Mm,n{C). (1.4) 



As announced in the introduction, our goal is to establish a LDP for {fJ.^)N for which 
the rate function involves a functional associated to a vector equilibrium problem, 
which itself describes the asymptotic distribution of the zeros of MOPs. Let us now 
explains our intuition. 

As it is classical for many random matrix models, one can embed such a non- 
centered Wishart matrix in a matrix valued stochastic process. Its squared singular 
values then induce (up to a time change of variable) a process of non-intersecting 
squared Bessel paths conditioned to start at a > and end at the origin. Kuijlaars, 
Martmez-Finkelshtein and Wielonsky studied the particle system of a fixed-time 
marginal and established a determinantal point process structure related to MOPs 
[38], a so-called MOP ensemble [36]. Moreover, the limiting zero distribution of the 
MOPs involved in this particle system has been characterized by a vector equilibrium 
problem in [39]. Combining these results, it is likely (see also [38, Appendix]) that 
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the spectral measure //^ converges almost surely as — t- oo to a limiting distribution 
^* which is the first component of the unique minimizer (/x* ,1^*) of the functional 

log^ — - — :dfi{x)dfi{y) - 1 1 log- — - — -dfi{x)du{y) 



F — y| J J \x — y 

+ 



when the vector of measures (/i, i^) runs over A^i(M-i-) x A^^y2 (^-)' where we in- 
troduced the set of constrained measures 

Ml/2 {^-) = ji^ G Mii2{^-) : du{x) < dx, ^(x) < ^ (i_g) 

Here we use the notation 

]R_ = (-oo,0], M+ = [0,+oo). 

Note that the functional (1.5) is actually not well-defined for all (/i, i^) G 7Wi(M+) x 
■M1/2 (^-)' since the logarithmic energy (1.1) can take the values +oo and — oo as 
well. We actually describe later an appropriate way to extend (1.5) to the whole 
set A^i(M+) X ^1/2 (^-)' which is possible because it lies in the class of weakly 
admissible vector equilibrium problems introduced by the authors in [32]. 

1.3 Statement of the result 

The aim of this work is to show that such a functional (1.5), once properly extended, 
is involved as a rate function governing a LDP for the spectral measures {fi^)]\f. 
More precisely, our main result is the following. 

Theorem 1.1. The sequence of measures (/i^)Ar satisfies a LDP on A^i(M-|_) in the 
scale N"^ with good rate function 

inf Ji-^u)— min J", 



|/iGA^i(M+): inf J(/i, i^) < 7) 



where J is a well-defined extension of (1.5) introduced in Section 3.1. Namely, 

(a) The level set 

[u^MA^A^): inf 

1-/2 (K_) 

is compact for any 7 E R. 

(b) J admits a unique minimizer on A^i(M-i-) x Ml/2^ 

(c) For any closed set T C 7Wi(IR+), 



limsup-^logPjvf/x^ G -F) < - inf [j{ii,v) - J{ii\v*)\. 
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(d) For any open set O C Mi{^+), 



lim inf 



1 



logPjvf^'^eo) >- inf {j{ii,u)-J{^\v*)\. 



As a direct consequence of Theorem 1.1 (b), (c) and the Borel-Cantelli Lemma, 
we obtain the almost sure convergence of towards //* in the weak topology of 
A^i(M+). Namely, if P denotes the measure induced by the product probability 
space (A^M,Af(C),PAr), we have 

Corollary 1.2. 

¥ (^11^ converges as N ^ oo to fi* in the weak topology o/A^i(M+)^ = 1. 
1.4 Generalizations and variations 

We now describe a few other particle systems for which one can use the same ap- 
proach as presented in this work to obtain a similar LDP statement. 



1.4.1 More general potentials 

The following generalization of the density distribution (1.2) has been introduced 
by Desrosiers and Forrester in [20] 

' (1.7) 



Z 



-NTr(y{X*X)-Re{X''A)) 



M,N 



where V : M4- — t- R is a continuous function which is extended to Hermitian matrices 
by functional calculus. Indeed, by choosing V{x) = x we recover the non-centered 
Wishart matrix model. Now, if we take A as in (1.4) and we assume that V satisfies 
the growth condition 

limmf — ^-^1 — — ^ > 1, 

x-y+oo 21og(x) 

then one can follow the methods developed in this work without substantial change 
to show an analogue of Theorem 1.1 and Corollary 1.2, where J is replaced by a well- 
defined extension in the sense of [32] (see also Section 3.1) over A^i(M+) x M.y^{W-) 
of the functional 

log - — - — -dfi{x)dfi{y) - / / log - — - — rdfi{x)du{y) 



\^ — y\ J J \x — y\ 

+ 



JJ log I — - — -^di^{x)di^{y) + j {y{x) -2y/lix^diJL{x). 
1.4.2 Rescaling the parameter a 

Observe that in our setting we have M/N — )• 1 as — )• 00. A natural question 
would be to investigate the case where one performs the rescaling a 1— )• aN , so 
that M/N — )• 1 + aasA^— T-oo with a > 0. It turns out that the approach we 
develop below is still well-suited for this case, but requires more involved asymptotic 



5 



estimates for Bessel functions and its zeros than the ones we use in this paper. These 
asymptotic estimates are actuahy provided by [1, (9.7.7)] and [1, (9.5.22)], and they 
would lead to statements similar to Theorem 1.1 and Corollary 1.2, where J is 
replaced by a well-defined extension of 

log^ — — rdfi{x)dfi{y) - 1 1 log- — - — :dfi{x)du{y) + [[ log- — - — -dv{x)dv{y) 
\x-y\ J J \x-y\ J J \x-y\ 



+ J yx — a log(2;) — \/ Aax + + a log ya + y/ 4ax + a'^j j d^{x) 
where the space -Ml 12 {^-) is now defined by 

dv Y^4a|x| — a 



A^^/2(M_) = <^zvGA^i/2((-oo,-f^]) : dv{x)^dx, -£{x) < 



2 



dx 27r\x\ I 

This is also the functional obtained in [39] when describing the limiting zero dis- 
tribution of the associated MOPs. Nevertheless, in this setting the proof becomes 
more technical, and we chose to restrict ourselves to the non-rescaled model for the 
sake of clarity. 

1.4.3 Non-intersecting Bessel paths vi^ith one positive starting and end- 
ing point 

In [18], Delvaux, Kuijlaars, Roman and Zhang investigated a system of N non- 
intersecting squared Bessel paths conditioned to start from a > at time t = and 
to end at 6 > when t = 1. It is actually not known if such model is related to a 
random matrix ensemble. We note that at fixed time < t < 1, it is easy to express 
the particle distribution as the marginal distribution of a Coulomb gas involving 
three different type particles by combining [18, Section 2.5] with the computations 
we present in Section 2. As a consequence, if we introduce the functional J to be 
the well-defined extension of 

log-: — - — 7dfi{x)dfi{y) - [flog-, — - — -dii(x)du{y) (11 
\x-y\ J J 

log- -dfi{x)dT]{y) + 1 1 log- -di'{x)di'{y) 




+ 1 1 TZ-—^]^dr]{x)d'q{y) + J l^J^]-—^ " " JZJ j ^^^^^ 

in the sense of [32] (see also Section 3.1), where /x G 7Wi(M+), and i^,r] £ 
satisfy 



^ix)<^\x\-y^ ^(x)< ^ 



1x1-^2 



dx nt ' dx 7r(l — t) 

then a LDP similar to the one of the non-centered Wishart matrix holds where the 
rate function is given by ^7 — min^ after taking the infimum over all constrained 
measures and ij. Indeed, there is no interaction between the particles associated 
to and rj, and then both u and rj interact with /i exactly in the same way that 
v interacts with /i in the present work, so that a LDP can be established with no 
extra work from the ingredients of the proof we present below. 
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1.5 Open problems 

There are other matrix models for which it is estabhshed that the Umiting mean 
spectral distribution is characterized in terms of the solution of a vector equilibrium 
problem, thanks to their connection with MOPs and a Riemann-Hilbert asymptotic 
analysis. Examples can be found in the Hermitian matrix model with an external 
source [10] and the two-matrix model [25, 23]. Nevertheless, it is not clear to the 
authors how to strengthen such convergence results to get LDPs. 

Another question of interest would be to see if the rate function introduced 
by Cabanal-Duvillard and Guionnet in [14] reduces for such matrix models to the 
functional of a vector equilibrium problem. 

1.6 Strategy of the proof 

In Section 2, we show that the joint eigenvalue distribution of the non-centered 
Wishart matrix is the marginal distribution of a 2D Coulomb gas with two type 
particles. The first type of particles are living on M+ and are exactly the eigenvalues 
of our matrix model. The second type of particles are abstract ones and live on a A^- 
dependent discrete subset of M_. They moreover attract the first type of particles, 
expressing the effect of the perturbation. This provides an insight as to why a 
functional like (1.5) should be involved as a rate function. To prove such statement, 
we first describe in Section 2.1 the eigenvalue distribution as a MOP ensemble, and 
then make use of the Nikishin structure satisfied by the weights associated to the 
polynomials in Section 2.2. 

In Section 3, we investigate the generalized particle system of the whole Coulomb 
gas for which we state a LDP, see Theorem 3.4. Theorem 1.1 then follows by 
contraction principle, as described by Corollary 3.5. We define the rate function in 
Section 3.1, which is a proper extension of (1.5), by following the approach of [32]. 
Prom the discrete character of the particles on M_ , a discussion provided in Section 
3.2 explains why the constraint set -A4^/2 (^-) naturally appears in the variational 
problem. 

In Section 4, we provide a proof of Theorem 3.4. The two main difficulties are 
the absence of confining potential acting on the particles living on M_, and the 
possible contact of the two different type of particles at the origin. Concerning the 
lack of confining potential, we follow an approach developed in [31] by one of the 
authors and perform a well-adapted compactification procedure. For the contact 
at the origin, we isolate the induced singularity and use the discrete character of 
the particles on M_ to control it, see the proof of Proposition 4.1 and particularly 
Lemma 4.4. 

Remark 1.3. From now, we will assume that N is even to simplify the notations and 
the presentation, but our proof easily adapts to the general case by replacing N/2 
by [A^/2] or [N/2\, and also 1/2 by [iV/2] /A^ or [A^/2j/iV, where it is necessary. 

2 A 2D Coulomb gas of two type particles 

In this section we show that the joint eigenvalue distribution is the marginal distri- 
bution of a Coulomb gas having two types of particles. Such a representation follows 
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from a particular type of MOP ensemble structure satisfied by the eigenvalues that 
we describe now. 



2.1 Multiple orthogonal polynomial ensemble 



We first show that the eigenvalues form a MOP ensemble in the sense of [37], a 
particular type of Borodin's biorthogonal ensemble [13]. For that, introduce the 
Vandermonde determinant 



An(x) 



det [a 



_niv 



J jj=i 



n 



(2.1) 



l<i<j<N 

For a > and a > 0, consider moreover the weight function 



-Nx 



X e 



where we introduced the modified Bessel function of the first kind 

^"^""^ "^^'r(A; + a + l) (2) 



(2.2) 



(2.3) 



We mention that these weights have been introduced and studied by Coussement 
and Van Assche [16, 17]. We now prove the following. 

Lemma 2.1. The joint probability density for the eigenvalues x = (xi, . . . ,xm) £ 
of X*X, when X is drawn according to (1.2), is a {N/2, N/2)-M0P ensemble 
with weights Wa,N CLnd Wa+i,N, that is given by 



—AN{x)det 



\x] '^a,N{Xj)^ 

\x) Wa+l,N{Xj)\ 

\ •' ) i,i=l . 



(2.4) 



where is a new normalization constant. 

Proof. We perform a singular value decomposition of X, that is we write X 
J7iXdiag U2 for unitaries Ui £ ZYm(C) and U2 G Un{C) with 



lag 



/XN 



Or 



and note that 

Tr {{X - A)*{X - A)) = Tr{X*X) - Tr{XA* + AX*) + Na 



(2.5) 



N 



J]^;, - Tr (Xdiagt/2vl*C/r + (Xdiagf/2A*[/i*)*) + Na. 



i=l 
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By integrating over the unitary groups, it follows from the Weyl integration formula 
[2, Section 4.1] and (2.5), that the probability density for the Xj's induced by Fj\f is 
given by 



fJl Jum 



NTr 



(^diag C/A* V* + (Xdiag V* ) * ) 



Um{C) JUn{C) 



dUdV, (2.6) 



where dU (resp. dV) stands for the Haar measure of WAr(C) (resp. Wm(C)) and Z^r 
is a new normalization constant. Note that one can assume the Xj's to be distinct 
since this holds almost surely. Consider 



B 



Or 



G MM,NiC) 



with a <bi < ••• < bi\i < a + 1. Then we have the following formula for the matrix 
integral [45, Section 3.2] 



f f ^NTr{x^,,,UB*V* + (Xm.sUB'V*r)^jj^y 

Jum{€) Jum{C) 



N 



det 



c„ n 



Ia{2N,/b~^ 



N 



AN{x)AN{b) 



, (2.7) 



where cn is a positive number which does not depend on x nor b. By continuity of 
the left-hand side of (2.7) in the fej's, we then obtain that (2.6) is proportional to 



Ajv(a;) 



lim • • • lim 



det 



and thus to 



Aj^{x) det 



db- 



i— 1 , 
I o 



a/2 



Lj2N^/bx 



b=a 



N 



(2. 



by I'Hopital Theorem. Finally, using for x > the relations [44, P79] 

^laix) = la+lix) + -laix), -^/q+i(x) = /^(x) - " ^ ^ Ia+l{x), 

dx X dx X 

it is easily shown inductively that the linear space spanned by the functions 



X I— )• 



-{x°/2J^(27V^/te)e-^"} 



b=a 



1 



, . . . , J. T 5 



matches with the one spanned by 



X I— )■ ^Wa,N{x), X I— >■ X* ^Woi+l,N{x) 



,N/2. 



This ends the proof of Lemma 2.1. 



□ 



9 



Remark 2.2. Although the parameter a associated to the matrix model is a non- 
negative integer, the distribution (2.4) still makes sense for non-negative real a. In 
fact, in the proofs we provide later, it will not matter whether a is an integer or not. 
Thus, if one considers the measures yJ^ (1-3) associated to Xj's drawn according to 
(2.4) with real a > 0, the LDP from Theorem 1.1 continues to hold. 



2.2 Nikishin system 

We now describe a property satisfied by the weights Wa^N and Wa+\,N-, a so-called 
Nikishin structure, and obtain as a consequence an exact Coulomb gas representation 
for the eigenvalues, see Proposition 2.4. The reader curious about Nikishin systems 
should have a look at [40] (where they are called MT systems). 

More precisely, it turns out that the ratio of the weights is (almost) the Cauchy 
transform of some measure, a fact which has already been observed [16, Theorem 
1]. We now make this result slightly more precise with an alternative simple proof. 
Consider the sequence 

< ^0,0 < Ja,l < jQ,2 < • • • 

of the positive zeros of the Bessel function of the first kind Jq,, a rotated version of 
Iq,, i.e 

J^(3;) = e^™/2/«(-ix), XGM+, (2.9) 
and introduce for each N the sequence of negative numbers 



2 



We then set for convenience 

Ajv = {afc,^ : > o} (2.11) 
and consider the associated normalized counting measure 

= ^ E '^(^)- (2.12) 

The weights Wa,N and Wa+i^^ then satisfy the following relation. 
Lemma 2.3. For all a > and a > 0, 



Wa,N VCL J X — U 

Proof. Up to a change of variable, this relation is nothing else than the Mittag-Leffler 
expansion 

which is itself provided by [27, P61] together with the relation (2.9). □ 
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Lemma 2.3 is in fact the key to express the eigenvalue density (2.4) as the 
marginal distribution of a two type particles Coulomb gas. Namely, if we introduce 
a Vandermonde-like product for (x, u) E 



N N/2 



N N/2 



^N,N/2{x, ^^) = n n ~ = n n 

i=l j=l 1=1 j=l 



(2.13) 



then the following Proposition holds. 

Proposition 2.4. The probability density (2.4) admits the following representation 



1 



a2 (^\a2 N N 2 

^N{^)^N/2y^) YJ r ^TTl u ( 

[[Wa,N{Xi) [[ \Ui\daN{ 



Ui 



1=1 



^N,N/2(X,U) 

where Zj^ is a new normalization constant. 

The proof we present now is inspired from the proof of [16, Theorem 2]. 
Proof. Recall that the density (2.4) is proportional to 

-1 N/2,N 



A AT (a;) det 



< X) Wa,N{Xj) \ 

J i,j=l 

<x) Wa+l,N{xj)> 

^ J t,j=l . 



(2.14) 



We first perform the factorization 

N/2,N 



det 



r _i ■\N/2,h 

\x] Wa,N{xj)\ 

) i,j=l 



1 N/2,N 
ii,j=l . 



det 



I ^ Wa,N ii,j=l 



N 



i=l 



and then use Lemma 2.3 to obtain 

1 N/2,N 



(2.15) 



det 



Provided with the identity 



'N/2 



N/2 



det 



^_^,N/2,N 



^ J i,j=l 



J 



X J ^^7, 



N/2,N 



N/2 

Y[daN{ui). (2.16) 

i=l 



X 



i-1 



X — U X — U 



X^^i-k+l 



k=0 
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the multilinearity of the determinant gives 



det 



{^-};r 



N/2,N 



■ det 



N/2,N ■ 

i,j=l 
\ N/2,N 



X I tl/rl 



:det 



-1 N/2,N 



N/2,N 



N/2 
1=1 



(2.17) 



Now, the well-known identity for mixed Cauchy-Vandermonde determinant, see e.g. 
[16, Lemma 3], yields 



det 



1 



N/2,N 



^N,N/2{X,U) 



(2.18) 



where the sign only depends on N. Combining (2.14)-(2.18), we obtain that (2.4) 
is proportional to 



N/2 



N/2 A 



N,N/2 



{x,u) 



Y[ Wa,N{Xi) JJ u\ d(JN{Ui 



(2.19) 



By summing the integrand of (2.19) over all possible permutations of the u[s and 
using the definition (2.1) of the Vandermonde determinant, we obtain that (2.19) is 
proportional to 



Al{x)Al^,{u) ^ 



N/2 



1^/2 An^N/2{x,u) f-J^ 

Since for every {x,u) G M.'^ x A^^^^ the quantity 



JJu;Q,,Ar(xi) Y[ \ui\ daN{ui). 



1=1 



Al^^{x)Ali^{u) 



N 



^N,N/2 



N/2^ 

'x,u] 



N/2 



1=1 



1=1 



is non-negative (and not identically zero), the new normalization constant Zjv has 
to be positive. The proof of Proposition 2.4 is therefore complete. □ 

In the next section, we perform a large deviations investigation for the whole 
Coulomb gas system. 
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3 A LDP for the generalized particle system 



On the basis of the preceding analysis, we investigate in this section the probabihty 
distribution on x M^^^ 

where, with u)c[,7v defined in (2.2), we introduced for convenience 

Vn{x) = -^logWa,N{x), X e M+. (3.2) 

The measure has been defined in (2.12), and Z^r is a normahzation constant. 
Consider the empirical measure for the second type particles 

N/2 



N 

i=l 



where the tij's are distributed according to (3.1) and note that the random vector 
of measures {^^ ,u^) takes values in A^i(M+) x M.ii2{R-), that we equip with 
the product topology. Our aim is to establish a LDP for ((^^,z^^))^, from which 
follows a LDP for {pl^)n by contraction principle. 

We first introduce the rate function in Section 3.1. Then, because of the discrete 
character of the second type particles, we introduce in Section 3.2 a convenient closed 
subspace of A^i/2(K-) where the u^^s actually live. Finally, we state the LDP for 
{lj,^,i'^)]\f in Section 3.3, see Theorem 3.4, and provide a proof for Theorem 1.1. 
The proof of Theorem 3.4 is deferred to Section 4. 



3.1 The rate function 

Our first task is to extend properly the definition of the functional (1.5) to A4i (M+) x 
1/2 (IK-)- A general method to do so has been presented in [32] and is based on a 
compactification procedure that we present for our particular case now. 



3.1.1 Compactification procedure 

Let S be the circle of centered in (0,1/2) of radius 1/2 and T : M ^ S the 
associated inverse stereographic projection, namely the map defined by 

\l + X'^ 1 + X"^ J 

It is known that T is an homeomorphism from M onto S \ {(0, 1)}, so that (S, T) 
is a one point compactification of M. For a measure on R, we denote by T^,// its 
push-forward by T, that is the measure on S characterized by 

[ f{x)dnfi{x) = [ f{T{x))d^x{x) (3.4) 
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for every Borel function / on §. We denote the two half-circles 

S± = {t(x) : X G M±} U {(0, 1)}. (3.5) 

Since T is an homeomorphism from IR+ (resp. ]R_) to \ {(0, 1)} (resp. S_ \ 
{(0, 1)}), it is easy to see (cf. [31, Lemma 2.1]) that T^, is a homeomorphism from 
A^i(R+) to 

{^JieMl{^+)■. M{(o,i)})=o}, 
and also from M.ij2{^~) to 

{/xGXi/2(S-): M{(0,1)})=0}. 

Equipped with such a transformation T^,, we are now able to provide a proper 
definition for the functional (1.5). 

3.1.2 Definition of the rate function 

Introduce the lower semi-continuous function V : S+ — s- K U {-I-cxd} by 

V(T(x)) =x-2Vax-^log(l + x2), x e M+, (3.6) 

and 

V((0, 1)) = liminf V(T(x)) = +oo. (3.7) 



We naturally extend the definition of the logarithmic energy (1.1) to measures on 
§ C (where | • | stands for the Euclidean norm) and define the functional J on 

A^i(M+) x7Wi/2(lR-) by 

J{li,v)= 11 log ^ dT^fi{z)dT^fi{w) - j j log ^ dT^^{z)dT^v{^) 

^ JJ^°^ J^^^dni^iOdTMC) + J V{z)dT,fi{z) (3.8) 

when both T^fi and T^u have finite logarithmic energy, and set J{^^v) = -t-oo 
otherwise. 

This definition is motivated by the following observation : from the metric rela- 
tion [4, Lemma 3.4.2], 

\T{x)-T{y)\= rr^t 2 ' ^^-^^ 

Vl + x^Vl + 

we obtain with (3.4) for any Borel measures /x, on M the relation 

[[ log -j 7dT^ii{z)dT^.v(w) = ((log-. Tdii(x)dv{y) 

JJ \z-w\ JJ \x-y\ 

+ ^^(E) j \og{l + x'^)dv{x) (3.10) 
+ ^i/(M) j log(l + x')(i^(x). 
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as soon as one assumes all these quantities to be finite. In this case, we obtain that 
J{^^v) matches with (1.5). 

Then, the following Proposition is a consequence of [32, Theorem 2.6]. 

Proposition 3.1. 

(a) The level set 

G Mi{R+) X 7Wi/2(M_) : J{f,,u) < 7} 

is compact for all 'j GM. 

(b) J is strictly convex on the set where it is finite. 

Because of the discrete character of the Uj's, we need to discuss now several 
constraint issues. 



3.2 Discreteness and constraint 

In this section we use the discrete character of the particles on M_ to build a closed 
subset £{R^) of A^i/2(M_) such that G £{R-) for ah N. This wih provide 
an explanation on why the measures on IR„ are restricted to the set Ai'^^^i^-) 
the minimization problem (1.5), and moreover will be of important use to control 
the possible contact at the origin of the different type particles during the proof of 
Theorem 3.4, see Lemma 4.4. 

We say that a measure v E A^i/2(IR-) is constrained by a Borel measure A on 
R_, that we note v < X, ii the signed measure A — is in a fact a (positive) measure. 
Introduce the set of constrained measures 



7W^/2(I^-) = { G Mi/2{R-) : < a} (3.11) 



and note it is closed. Indeed, if (^'Af)Af is a sequence in ^Ay20^-) with weak limit 
v, then (A — i^Ar)7v converges in the vague topology (i.e the topology coming from 
duality with the Banach space of compactly supported continuous functions on M) 
towards X — v, which is hence not signed. 

Since the random variables Uj's take values in A^r, see (2.11), we have almost 
surely 

1 ^^Z' 1 

i=l mGAjv 

and thus almost surely G -^1/2 (^-) ^'^^ ™y Consider the measure a on M_ 
having for density 

^ = ^1x1-^2. (3.12) 

ax TT 

and note that the Radon-Nikodym theorem yields that the definition of A4y2(M-) 
presented in this section matches with (1.6). It is in fact the limiting distribution of 
the constraints (Tat. 

Lemma 3.2. The sequence {(Ti\i)m converges towards a in the vague topology. 
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Proof. Since for any 6 < we clearly have limiv->oo ^Af ({&}) = (^{{t>}) = 0, it is 
enough to show that for all 6 < 0, liniTv-s-oo (^N{[b, 0]) = a{[b, 0]), that is 

lim -^^\k : Uk^N > ^| = — / 

Af->oo N I J Jb 

By change of variables, it is equivalent to prove that for all 6 > 



dx = —. 



AT-i-oo iV L iV J TT Jo 

Fix e > and let k{N) be the integer part of ^^'^^^ . The McMahon asymptotic 
formula [1, formula 9.5.12] yields 

lim = TT, (3.13) 

k^oo k 

and thus 

lim '^ = b + s. 
As a consequence, we obtain the upper bound 

limsuT) —ii k : < b\ < lim = — ^ 



Similarly, changing e by — e in the definition of k{N) yields the lower bound 

liminf — tjlfc : < b\ > 

and Lemma 3.2 follows by letting e — )• 0. □ 

We now introduce the subset i£'(R_) of A4i/20^-) of the measures which are 
either constrained by a^, for some A^, or by a. Namely, 

oo 

£{R.)= U Mlf,{R-)[jM^y,{R.). (3.14) 

N=l 

By construction u'^ E iS(M_), for any A^, and moreover 
Lemma 3.3. (^(M^) is a dosed subset of A4i/2{R~)- 

Proof. Let be a sequence in i5(M_) with weak limit i/, and let us show that 

1/ G £{R-). Since the sets Ml^^ C^-) -^1/2 (^-) closed for all A^, one may 
assume that i^j < cr]\f., with liuij^oo = +oo. One then obtains by Lemma 3.2 
that u < a, and thus u e iS(M_). □ 

Concerning the measure on S_, see (3.5), we similarly set 

oo 

£{8.) = [j M^'/^-{8^)(jM^J^{S.), (3.15) 

N=l 

SO that T^v^ G "^(S-) for any A^. Moreover, note that since ^^({(0, 1)}) = for any 
V E £^(§_), it follows that T* is an homeomorphism from (^(M-) to iS(S_), and <£"(§-) 
is seen to be a closed subset of 7Wi/2(S-) from Lemma 3.3. 
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3.3 LDP for the generalized particle system 

We are now in a position to state the LDP for {^j,^ ,v^)n ■ Let us precise that we 
equip f (M_) with the topology induced by A^i/2(IR-) and 7Wi(M+) x £'(M_) carries 
the product one. Then the fohowing LDP holds. 

Theorem 3.4. The sequence {fi^,u^)N satisfies a LDP on x £r(M_) in 

the scale N"^ with good rate function J — min JT". More •precisely, 

(a) The level set 

G X f (IR_) : J[ii,v) < 7} 

is compact for any 7 G M. 

(b) J admits a unique minimizer {fi*,^*) on Mi(K-^-) x 

(c) For any closed set T C A^i(M+) x <S(M_), 

limsup-^logFjvfl;"'^,/^'^) G^) < - inf | z.) - 

(d) For any open sei O C x £:(M_), 

liminf-^logP^((/.^,i^^)GO) >- \ni {J{^l,y)-J{,,*y)]. 

A direct consequence of Theorem 3.4 is Theorem 1.1. 
Corollary 3.5. Theorem 1.1 holds true. 

Proof. Theorem 1.1 follows by contraction principle (see [19, Theorem 4.2.1]) along 
the projection A^i(R+) x f(]R_) — ^ A^i(M+) and the fact that J{ii,u) = +00 as 
soon as 1/ G \ Mli^i^..). □ 

4 Proof of Theorem 3.4 

We first observe that Theorem 3.4 (a), (b) easily follow from Proposition 3.1. 

Proof of Theorem 3.4 (a), (b). Since £^(M_) is a closed subset of M.i/2{^~) (see 
Lemma 3.3), Theorem 3.4 (a) follows from Proposition 3.1 (a). The existence of a 
minimizer for J' on A4i/2(M-) ^ <f (IK__) is a consequence of Theorem 3.4 (a). Since 
the set My20^-) is convex, and i7(/i, = +00 as soon as G <?(M_) \ My2(^-): 
the minimizer is unique by Proposition 3.1 (b). □ 

Concerning the proof of Theorem 3.4 (c), (d), it is usually pretty standard to 
establish LDP upper and lower bounds by proving a weak LDP and an exponential 
tightness property (see [19] for a general presentation on LDPs). However, because 
of the lack of confining potential acting on the particles on M_ , it is not clear to the 
authors how to prove directly that the sequence {ix^ is exponentially tight. 
Instead, we follow a different strategy developed in [31] : we first prove in Section 
4.1 a weak LDP upper bound for {T^fi'^ ,T^v^)j\i, the push-forward of {^j,^ ,v^)n 
by the inverse stereographic projection T. We then establish a LDP lower bound 
for {ji^ ,v^)i\i in Section 4.2, and show in Section 4.3 that it is enough to obtain 
Theorem 3.4 (c), (d). 
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4.1 A weak LDP upper bound for [T^ii^ ,T^v^)n 
Consider the functional J on A^i(S+) x defined by 

J{^J',v) = jj log 1^ ]_ ^^ dn{z)dfi{w) - jj log ^ ^^ dii{z)dv{i) 



^ jj^''^ ^^\Mi)d^{0 + j V{z)di^{z) (4.1) 



if both /i and v have finite logarithmic energy, and set J(/i, v) = +oo otherwise. We 
recall that V has been introduced in (3.6)-(3.7) and £"(§-) in (3.15). Note that, 
with J defined in (3.8), the following relation holds 

Jill, v) = Jinn, T^v), (fi, i^) G Mi{R+) X S(R_). (4.2) 

Now, choose a metric compatible with the topology of 7Wi(S+) x£{S-) and write 
Bsin, v) for the open ball of radius 5 > centered at (^, v). The aim of this section 
is to establish the following weak LDP upper bound for [T^u^ ^Ti^v^). 

Proposition 4.1. For any {im^v) £ Mi(S+) x f (S-) 

limsuplimsup-^log(zArPArf (r*^^,r*z^^) e BsifJ.,i'))\ < -J{fi,u). (4.3) 

Concerning the proof, we first describe in Section 4.1.1 the induced distribution 
for the particles {T{xi)).^-^ and (T{ui)) on X %^J\ Then, we show (4.3) in 
Section 4.1.2, where the main difficulty is to control the singularity created by the 
fact that the different type particles may meet at the origin when N ^ oo. To do 
so, we will use a few technical lemmas, for which the proofs are deferred to Section 
4.1.3 for convenience. 

4.1.1 The induced distribution for the particles on S 

Introduce the random variables on S 

Zi = T{xi), i = l,...,N, ^^ = T{ui), i = l,...,iV/2, (4.4) 
where the Xj's and the UiS are distributed according to (3.1). Thus 

. N N/2 
1=1 i=l 

We set the measures A = r*(l]R_^(x)drc) on S+ and r/7v = T^ctn on S_, with 
introduced in (2.12). From V/v introduced in (3.2), we also construct the lower 
semi-continuous function Vat : S+ — )• M U {+00} by 

VN{Tix)) =VN{x)-^logil + x^), XGM+, (4.6) 

and 

Viv((0,l)) =liminfViv(r(x)) =+00, (4.7) 
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where the latter equaUty fohows from the asymptotic behavior [1, formula 9.7.1] 



Ia{x) 



V2' 



TTX 



1 + O ^) j as X — )• +00. 



(4.8) 



Then the following holds. 

Lemma 4.2. The joint distribution of {z, ^) = (zi, . . . , zat, ^i, . . . , C7V/2) ^■^ given by 



1 



^Af,Ar/2(2,^) 



N 



N/2 



■t=i 



'^dX{z,)ll\C^\Vl^^dm{C^) 



1=1 



where Zj^ has been introduced in (3.1). 

Proof. From the metric relation (3.9) we obtain 

N 



A^(^) = |A^(r(^))in(i + ^?) 

N/2 



2\N-1 



A%/,{u) = \A%/,{T{u))\l[{l + u^) 



2 \ N/2-1 



i=l 



N 



N/2 

^N,N/2{x,u) = \Ar,,^/,{T{x),T{u))\ll{l + x^f^'l[{l + u^f/\ 



1=1 



i=l 



Thus, with Vat defined in (4.6), this yields 



N 



N/2 



Ui 



^N,N/2ix,u] .^^ .^^ 

A%{T{x))A%^,{T{u)) 



An^n/2{T{x),T{u)) 

We moreover obtain from (3.9) the identities 
1 . .,0 |x| 



11 1 + ^2 111, 



i=l 



i=l 



(4.9) 



1 + x2 

and then from (4.9) 



1 - |r(x) 



1 + 3;2 



|r(x)|Vl-|r(x)|2, xEM, (4.10) 



^%ix)A%/2i^) 



N 



N/2 



^N,N/2{X,U) 

A%{T{x))A%^,{T{u)) 



Ui 



i=l 



An,n/2{T{x),T{u)) 



N 



N/2 



n (i-ir(x.)ne-^^-(^(-))^ n m-oi v^r^"^ 



i=l 



i=l 



(4.11) 



Lemma 4.2 then follows from (4.11) by performing the change of variables Zi = T{xi) 
for i = 1, . . . , iV and = T{ui) for i = 1, . . . , N/2. □ 
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4.1.2 Core of the proof for Proposition 4.1 

Provided with Lemma 4.2, we now establish Proposition 4.1, up to the proofs of few 
lemmas which are deferred to the next section. 

Proof of Proposition 4-1- We obtain from (4.5) and Lemma 4.2 



N 

X 

i=l 



N 



-NVNizi) 



i=l 



^N,N/2{Z,^) 

N/2 

Ylil - \zi\^)dX{zi) H mV^-\^^?drJN{^^). (4.12) 



i=l 



We write 



exp 



N 



-Qg-iVViv{^,) 



i=l 



N N/2 \^ 
+ {2V N{Z^)+ log \zi-Cj\) I 

i=l 7=1 J / 



exp - 



log- 



N t 



\z — w\ 



(4.13) 



1 



+ / / log ^^dT,z.^(OdT,r.^(C) 



2VAr(z) + log \z - ei )dT,fi^{z)dTy' (e) 



Note that, since T^iu.^ (E> T*/i^{(2;, to) G S+ x S+ : z = = 1/N almost surely, for 
any M > we have almost surely 



1 



log ■ 

z=iw \Z-W\ 



dT,fi^{z)dT,fi''{w) 



N ( 



> I I min ( log ^ MyT,n^{z)dT,n^{w) - ^ 



(4.14) 



and similarly 



iog-l-dr,i.^(e)dr,z.^(C) 



> / / min I log 



^l^,M)dr.i.^(Odr,^^(C) - ^. 



(4.15) 
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To make the control of the singularity at the origin easier, we write for any M > 
; 2Vn{z) + log \z - i\^dT,ii^{z)dT,v^{^) 
2Vn{z) + log \z-i\- log \i\]dT,^l^{z)dT,u^{i) 



+ j \og\i\d%u''{0 
> JJ mm (2VNiz)+log\z- (\-log\^\,MyT,fi^iz)dT,iy^i() (4.16) 
+ J log|C|dT,z.^(C). 

Note that the latter step makes sense since T^v^ can not have a mass point at (0, 0). 
Such a decomposition is motivated by the following lemma. 

Lemma 4.3. For any G NU {oo}, the map 

{z, ^ 2Vn{z) + log |z - - log ICI (4.17) 

is bounded from below on x S_, where we denote Voo = V. 

Now, if we introduce for any M > and (/u, u) G A^i(S+) x <?(§-) 

J^'if^,iy) mm[\og—^ My ^x{x)dii{y) (4.18) 



+ 



j j mm (2Vn{z) + log |z - C| - log m) d^i{z)du{i) 
+ 11 min(log^l^,M)dKO^^KC)+ I log|e|a!Ke), 



we obtain from (4.12)-(4.16) that 

ZNFN(iT,fi^,T,i^^)eBs{fi,i^)) <CNexp\ -N^ inf J^}, (4.19) 



where we set 

„ AT N/2 

= e3^^^/2 / - |z,|2)dA(z.) n UV^^driMii-d. 

Note that by construction J^^ is bounded from above, but may take the value — oo 
for some (//, z^) G Mi{E>+) x Mi/2{S-). Our choice to restrict A^i/2(§-) to <?(§-) 
is motivated by the following key lemma, which yields in particular that Jj(Jf is well 
defined and has each of its components bounded on A^i(S+) x <?(§_). 

Lemma 4.4. The functional 

'log\^\di^{0 
is continuous, and thus bounded, on <?(§-) 
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We observe that 
Lemma 4.5. 



lim sup log Cn < 0. (4.20) 



As a consequence, we obtain from (4.19) 

limsup-^log|zArPf (T,/x^,r,i^^) Ge5(/i,z^)) j < -liminf inf J^. (4.21) 
Now, introduce for any M > and (//, v) E A^i(S+) x <£"(§-) 

J^\^^v) = //minflog , ^ . ,M)dn{x)dn{y) (4.22) 



,z — 

+ ^y"min (2V(z) + log|z-e| -log|e|,M)d/i(z)dz^(0 

+ min(log^^,Af)(ii/(OMC) + y log|C|a!Ke) 

since the fohowing holds 
Lemma 4.6. 

lim inf inf J;Jf > inf . (4.23) 
It thus follows from (4.21) that 

limsnp ^log\ ZnF( {nfi^,ni^^) £Bsif^,M <- inf J*^. (4.24) 

Note that for any M > 0, the function 

(z, w) I— >• min ( log ■] r, M 

V — 

is continuous on S x S, so that the functional 

/i I— J J min ^ log 1 -, M^dfx{z)dfi{w) 

is continuous on A^i(S+), as well on Lemma 4.3 moreover yields for any 

M > the continuity of 

ifi,!^)^ jj min {2V{z) + \og\z - i\ - \og\i\,M)diJi{z)du{i). 

Thus, this shows with Lemma 4.4 that J^^ defined in (4.22) is continuous on 
A^i(S+) X <?(M_), and we obtain by letting (5 ^ in (4.24) that 

limsuplimsup-^loglZjvPf (T,/x^,r,i/^) e Bs{^i,v)] \ < -J^\n,u). (4.25) 
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Letting M — )• +00 in (4.25), the monotone convergence theorem yields 
hmsupUmsup-^log jzTvPf (T^fi^ ,T^u^) G Bs{lJ-,i^)]\ 

< - J 1 1 log ^ dfi{z)dfi{w) (4.26) 



+ JJ [2V{z) + log\z - C\ - log\C\)dfi{z)du{0 

+ 11 iog^^t^KOMC) + 1 loglcic^Koj- 

Finahy, in order to obtain Proposition 4.1 from (4.26), it is sufficient to show that, 
with J defined in (4.1), 

J(^, iy) = [[ log -^dfi{z)dfi{w) (4.27) 
././ \z-w\ 



+ JJ [2V{z)+log\z-C\-log\^\jdfi{z)du{0 

+ 11 ^ogj^^du{Odiy{C) + I log\C\du{C) 

for all (/U, I') £ 7Wi(S-|-) x £^(S_). Note that if or has infinite logarithmic energy, 
then Lemmas 4.3 and 4.4 yield that the right-hand side of (4.27) is +00. If both 
H and have finite logarithmic energy, then it is known (see e.g. [32, Section 3.1]) 
that 

/ [ log < +°°' 

and thus since V is bounded from below 

II (2V(z) + log|z-e|-log|e|)(iM-2)M0 + l^og\C\dHO 
= I V{z)di^{z) - // log ^l^d^(z)di.(0, 

which proves (4.27). The proof of Proposition 4.1 is therefore complete, up to the 
proofs of the lemmas. 

□ 

4.1.3 Proofs of Lemmas 4.3, 4.4, 4.5, and 4.6 

Proof of Lemma 4-3- We have the inequality 



i^-ei > leiVi- kP, ces^. (4.28) 

Indeed, (4.28) trivially holds if z = (0,1). Since for any z G S the Pythagorean 
theorem yields \z - (0, 1)| = y^l - |z|2, (4.28) moreover holds when = (0, 1). If 
none of z or ^ is (0, 1), then there exist x G M+ and u G M- such that |z — ^| = 
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\T{x) — T{u)\. Inequality (4.28) then follows from the metric relations (3.9), (4.10) 
and the inequality |x — ii| > \u\ when {x,u) G M+ x M_. 

As a consequence of the inequality (4.28), we obtain for any (z,^) € S+ x S_ 
and iV G N U {oo} 

2Vn{z) + log \z-^\- log lei > 2Vn{z) + \ log(l - \zf). (4.29) 

Now, from the the metric relations (4.10) we obtain 

inf f2Voo(z) + ^log(l-|z|2)) =2 inf fx - 2^^ - log(l + x^)) > -oo, (4.30) 

and similarly for any G N, 

inf f2V7v(^) + \ log(l - \z\^) \ = 2 inf (Vt^ix) - log(l + x^)] > -oo, (4.31) 

where the latter inequality follows from the definition (3.2) of V^r and the asymptotic 
behavior (4.8) of the Bessel function. Lemma 4.3 then follows from (4.29)-(4.31). □ 

Proof of Lemma 4-4- Since T* is an homeomorphism from (^(M^) to £{§-), we ob- 
tain with the metric relation (3.9) that for any v G <?(§-) 



\og\^\du{^) = [ log\T{u)\dn-'u{u) 



log ( , „ ) dT,-^u{u 



\u 



2 



/ log\u\dn-^u{u) + F{u), 

J\u\<l 



'I 

where F is a continuous function on £^(S_). Lemma 4.4 is thus equivalent to the 
continuity on i?(R_) of the functional 



/ log\u\di^{u), (4.32) 

J\u\<l 

which is itself equivalent to the uniformly integrability of u i— t- 1|„|<i log with 
respect to the measures of (?(]R_), namely to 

lim 



sup / |l|„|<i log |u| |(ii^(u) = 0. (4.33) 

ue£{R-) J\u\<e 
Since for any e > and any v G iS(M_) 

/ |l|n|<ilogKI \du{u) < I log^ \u\du{u), 

J\u\<e |log(e)| 7|„|<i 



^1 

it is enough to show that 



sup / log |npdz^('u) < +00 (4.34) 

ue£(M.-) J\u\<l 
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in order to obtain (4.33). By definition (3.14) of if (M_) we liave 

sup / log^ \u\du{u) (4.35) 

ue£{R-) J\u\<l 

< max < sup / log^ \u\daN{u) , / log^ \u\da{u) > . 

[ N J\u\<l -'|«|<1 J 

First, it follows from the definition (3.12) of a that 

/ log^ \u\da{u) = — [ x^/^log^{x)dx < +00. (4.36) 

J\u\<l ^ Jo 

Then, the definition (2.12) of gives 

It is a consequence of the McMahon expansion formula [1, formula 9.5.12] that 

lim (ja,k+l - ja,k) = TT, (4.38) 
fe— s-oo 

and this provides the existence of C > independent of satisfying 
1 , 2 f ja,k \ 



2 



< C [ log'^{x^)dx < +00. 
Jo 



(4.39) 



Indeed, the latter inequality follows by splitting the integration domain and from 
the fact that x i— )• log^(x^) is non-negative and decreasing on [0,1]. Combining 
(4.35)"(4.37) and (4.39) we obtain (4.34), which completes the proof of Lemma 4.4. 

□ 

Proof of Lemma 4-5. From the metric relations (4.10) we obtain 

„ N N/2 

CN = e'''/' / n(l-k.|'MA(z,)ni^^lA/l^ff^^^(^') 

= e3^/2Q (i-lzndA(z)) IClVT^dmio) 

/ r 1 \^ / r II \ 

-'"'^ X, r^^^- / r^L....(.) . (4.0) 
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Since the definition (2.12) of ajsi yields 

then Lemma 4.5 follows from (4.40) and the identity [44, Section 15.51] 

^ 1 1 

Proof of Lemma 4-6. We write 



□ 



liminf inf > inf J^+liminf inf (j^-J^) (4.41) 

and note that, from the definitions (4.18) and (4.22) of Jf^ and J*^ respectively, we 
have 

inf (4^-jA^)>i inf |min(2V7v(^)+log|2-C|-log|C|,M 

-minf2V(z) +log|z-C| -log|C|,M)|. (4.42) 



The inequality (4.29) and the fast growth of V(2;) and Vn{z) as z — )• (0, 1), which 
follows from the definitions (3.6)-(3.7), (4.6)-(4.7) and the asymptotic behavior 
(4.8), provide the existence of a neighborhood A^oo C §+ of (0, 1) such that for all N 



m\n[2V n{z) + \og\z - i\ - \og\i\,M^ (4.43) 
= min(2V(z)+log|z-e| -log|e|,M) =M, (z, € A^oo x S_. 
Next, we claim the existence of a subset A/q C S+ satisfying A/q U AAqo = S+ and 

min(^2Vjv(^) +log|2-C| -log|e|,M) (4.44) 
> minf2V(z) + log|z-C|-log|C|,M), (z, ^ G A/'o x S_, 



for any A*" sufficiently large, so that Lemma 4.6 would follow by combining (4.41)- 
(4.44). 

To show this it is enough to prove that for any L > there exists > such 
that for all N > Nl 

Vn{x) - X + 2y/ax > 0, xG[0,L], 

or equivalently (see the definitions (3.2) and (2.2)) 

y"/a(y)e"^ < (2iV^/^)", ye[0,2NV^]. (4.45) 

Indeed, if we choose A/q = T{[0,L]) with L large enough so that A/q U A/'oo = §+, 
then (4.44) would hold for any N > Nl as a consequence of (4.45). Given L > 0, 
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if a = then (4.45) holds because /o(0) = 1 and y i— t- Io{y)e ^ is decreasing on 
M+. If a > 0, it is then easy to see from the asymptotic behavior y"Ia{y)e~^ = 
(27r)-V2y«-i/2(i + O(y-i)) as y -)> +oo, provided by (4.8), that (4.45) is satisfied 
for any N large enough. This completes the proof of Lemma 4.6. 

□ 

We now provide a proof for the announced LDP lower bound. 

4.2 A LDP lower bound for (/i^,zy^)Ar 
The aim of this section is to establish the following. 
Proposition 4.7. For any open set O C 7Wi(R+) x £^(M_) 

hminf log jZjvPiv ( {f^"" ^i^"") ^ o)} > - inf J{fi, v). 

Proof. Note that it is sufficient to show that for all (//, i/) G O 

lim^nf ^ log [zMPN^ifi"", i^"") e O) } > -Jifi, y). (4.46) 

We first prove in two steps that (4.46) holds if and u satisfy the following : 
Assumption 4.8. 

(1) /i and V have compact support. 

(2) Supp(^) C M+ \ {0} and Supp(z^) C M_ \ {0}. 

(3) With a as in (3.12), there exists < e < 1 such that < (1 — e)a. 

(4) T*;U and T^v have finite logarithmic energy. 

We then extend in a last step (4.46) to all (/U, v) ^ O hy mean of an approximation 
procedure. This approach is similar to the strategy developed in [6, Section 3.2], see 
also [28, Section 3.4]. 

Step 1 (Discretization) Given (^, v) ^ O satisfying Assumption 4.8, our first 
step consists to build discrete approximations of (/x, v). To this aim, we note that 
and V have no atom as a consequence of Assumption 4.8 (d) and consider 

x\ = mm < a; G 



{xGR+i ^([0,x]) =^}, (4.47) 
min{x>xf^ ^([xf\x]) = 1}, i = l,...,7V-l, (4.48) 



and similarly 

yj"^) = min {y G M_ : z.((-oo, y]) = (4.49) 
yi^)=min{y>yf): z.([yf ), y]) = 1}, i = 1, . . . , Ar/2 - L (4.50) 
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Since ^ and v moreover have compact supports, the following weak convergence 
follows easily 

hm - )) = ^ and lim - b{yf^) = v. (4.51) 

i=l 1=1 

Because the Uj's are distributed on the discrete set Aat (2.11), we also set 

lij^^ = maxjn G Aat : u < yj^^j, i = 1, . . . , A^/2, (4.52) 
and moreover introduce 

N/2 
i=l 

We now show that, for any N large enough, the uf^^^s lie in the convex hull 
co(Supp(z^)) and the following interlacing property holds 

yr'<u^l<ySl ^ = !,■■■, N/2-1. (4.54) 
Indeed, with e as in Assumption 4.8 (3), (4.38) yields such that 

sup {ja,k+l - ja,k) < 7r(l + s) 
k>ks 

and, since ^ Supp(i^) by assumption, there exists N^^ such that 

sup iy[[ak+l,N, o-kM) =0' N > Ne. 
k<ke 

Thus, recalling the definition (2.10) of the a^^Ar's, we obtain for any N > Ni; 
supi^{[ak+i,N,ak,N]) = sup iy{[ak+i,N,ak,N]) 

k>0 k>ks 

< (1 - e) sup a[[ak+i,N, afc.Af]) 



1 

" k>k^ 



^1 - sup (ia,fc+l - ja,k) 

vriV ■ - ■ 



< fl-e 



1 



N 

The latter inequality implies that there exists an element of A^r in each {y\^\y\^\) 
provided is large enough, so that (4.54) follows from the definition (4.52) of the 
Uj's, and moreover that all the UiS are in co(Supp(i^)). 

Note that (4.54) yields v'^^^ < aN, and thus z^(^) G £{R-) for ah N. Moreover, 
by combining (4.54) with (4.51), we obtain the weak convergence of {i'^^^)n towards 
u. As the result of the discretization step, we have shown the existence of 5o > 
and A^o such that for all < 5 < 6q and N > Nq 

J^(^(xi),i^^^^) : xGM^, mLx\xi-xf^\<6^ CO. (4.55) 
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Step 2. (Lower bound) We now prove (4.46) when {fi, v) satisfies Assumption 
4.8. As a consequence of (4.55) we obtain for any {) < 5 < 5q 



r /S? ('r\ A2 (iiW^ N 

> A (X) A ) I W|-rT^-w.(..)^ 

— /r , (iv), ,.1 A,...,./^^ ,,(Af)^ 11' « '11 * ^ ^ 



i=l 



For a Borel measure A on M witli compact support, introduce its logarithmic poten- 
tial 

U^ix) = I log , ^ , d\{u) 
J \x — u\ 

which is continuous on M \ Supp(A) [42, Chapter 0] and note that 

N 

A^,^/2(a;,^(^)) =nexp{ -iVC/'^'"'(x,)}. (4.57) 

i=l 

We also set for x € M+ 

W7v(x) = y7v(x) (x), (4.58) 
W{x) =x-2y/^-U''{x) (4.59) 

and obtain from (4.56)^(4.59) 

z^p^((Ai^^)Go; 

> exp|-iV2 max |Tyiv(x) - VF(x)||A^/2(tx(^))|ao,ivr/^ (4.60) 

AT 



i=l 



By using the change of variables Xi^ Xi + x^^^ for i = 1, . . . ,N, and the fact that 
\xl Xj + Xi — Xj\ > max||Xj' — Xj \ , \Xi — Xj\j as soon as Xi > Xj and 
Xi > a;^- , we hnd 

N 



> / A^(a; + ^W)TTe-^^(^^+^'"')dx, 

Jm"" r=i 



> n (-f ^ - n i-iti - n ^-""^^'^"^ (^-^d 

i+Kj i=l i=l 

. N-1 N 



29 



Since the xj^^'s lie in the compact set co(Supp(/x)) and W is continuous there, we 
obtain 

hmhmsup max max \W {x + x[^^ ) — W {x[^^ )\ = (4.62) 

and also, using moreover (4.51), 

i ^ 
im — 

N 



1 ^ r 

lim — ^T^(xS^))= W(x)dfi{x). (4.63) 
i=i 



Using the change of variables ui = xi and lij+i = Xj+i — for i = 1, . . . , iV — 1, it 
follows 



N-l 



. dxiY\_ (xi+i - Xi)dxi+i 



> 



duiT\uidui = -j^ [ — ) . (4.64) 



1=2 

We thus obtain from (4.61)-(4.64) 

1 r ^ 
liminfliminf ^log / AI (x)Y\ e-^^^'^^^d 



N-l 



Next, we have 



^i+l<j i=l 

- I Wix)dfiix). 



lim max \Wn{x) - W{x)\ = 0. (4.66) 

N—^QO a;Gco(Supp(/^)) 



Indeed, the asymptotic behavior (4.8) yields the uniform convergence of Vi\i{x) to- 
wards x — 2-y/ax as — )• oo on every compact subset of M+ \ {0}, and in particular 
on co(Supp(^)). It is thus enough to show the uniform convergence ofU"' ' to 
on co(Supp(^)) to obtain (4.66). For any x G Supp(/i), the map y i-)- log|x — y\ is 
continuous and bounded on co(Supp(i^)), so that the pointwise convergence of JJ'^^^^ 
to on co(Supp(^)) follows from from the weak convergence of u^^^ to z^. Since 
for all N the map C/*^'^' is continuous and decreasing on the compact co(Supp(//)), 
and that is moreover continuous there, the pointwise convergence extends to the 
uniform convergence by Dini's theorem. 

We thus obtain from (4.65)-(4.66) by taking the limit N ^ oo and then (5 — )■ 
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in (4.60) that 

\i+i<i i=i / 

+ liminf ^ log [u^j^^ — uf^'^Y ~ j W{x)d^{x). 

i<j 

Now, note that because x i— )■ log{x) increases on M+ the definition (4.47)-(4.48) of 
the x-^'^'s yields 

± E log (.f - .rv + E log (42 - 4"') 

= 2 log {xf^\ - )) / / l.<, 

> 2 I / ^ log(y - (4.68) 

and then that 

2 lim / / log(w — x)dii(x)dii(y) 

log\x — y\dfi{x)dfj,{y). (4.69) 
The interlacing property (4.54) yields 



(N) (N) ^ (N) (N) f • , 1 / • 

Uj -Ui > Vj-i -Vi tor z + 1 < J, 



and thus 



7V/2-1 

5:iog(.f)-.s-))^> iog(c^-nS-))^+ ^og{yn-yrr- 

i<j 1=1 i+Kj 

Since 

min {u^l - uf^) > inf (a.,^ - a,^,^) 

l<i<N/2 ^ ^ ^ fe>0 ^ ' ' / 

- 2Sv^gi(^"''=+i-^''^''^)' 

we obtain from (4.38) and (4.70) 

31 



Moreover, because for 1 < i < N/2 — 1 

- ,f ) > 2| max(Supp(.))|V2(|,f )|V2 _ |yW|i/2 ) 



> 



^|max(Supp(.))|V2a(bf),y(f) 
V 

1/2 



^|max(Supp(.))|V2,([,^^-^y^^-/]) 
= ^^l™ax(Supp(i/))|^/^ 

we obtain from (4.71) 

'KiifME'°«^(^"'-4"')^>li.nmf-l log(.fi-!,f))l (472) 

i<j i+2<j 

Next, similarly than in (4.68)-(4.69), we obtain from the definition (4.49)-(4.50) of 
the yl^^^s that 

^tJ^^ w E (yj-l - yi^^)^ 



i+2<j 



= 2Hminf log — -uj^'') // lu<v du(u)dij(v) 

> 21iminf // log,(v — u)di'(u)di'(v) 



\og\x — y\dv{x)dv{y). (4.73) 
From (4.67)^(4.69) and (4.72)-(4.73) it follows 



> - I log ^ ^^ dfi{x)dfi{y) + j W{x)dii{x) (4.74) 

+ jj log 1^ \ y^ dv{x)dv{y)^. 

Since both V and are bounded and continuous functions on the compact Supp(/x), 
by (4.59) 



W{x)d^M{x) = j {x-2yja^^dn{x)- j U''{x)dn{x) 

= j — 2y/ax^diJL{x) — j j log 1 -^d^{x)dv{y), 



and thus 
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> - <( / / log - log (4.75) 



+ 



J J log j ^di'{x)di'{y) + J (^x — diJ,{x) 



Since by assumption the measures /j,, v have compact supports and T*;U, T^v have 
finite logarithmic energies, then /i, v also have finite logarithmic energies and clearly 



y log(l + x^)dix{x) < +00, j log(l + 



x'^)di'{x) < +00. 



Thus, one can use the relation (3.10) and obtain that the right-hand side of (4.75) 
equals J{^^v), see (3.8), which proves (4.46). 

Step 3. (Approximation) First note that (4.46) trivially holds as soon as 
J'(/i, z^) = +00. It is thus enough to show (4.46) when both T*/i and T*z^ have 
finite logarithmic energy, and one can moreover assume that u < a. For such (^, v), 
we now construct a sequence (^fc, v^jk of A4i(M+) x £^(R_) where each (/i^, v^) sat- 
isfies Assumption 4.8, such that we have the weak convergences 

lim yUfc = /u, lim = 

fc— >oo fc— ^-oo 

and which moreover satisfies 



lim J{f,k,i^k) = J{f^,i^). (4.76) 

This, combined with the two first steps of the proof, shows that (4.46) actually holds 
for all {fj,, u) G O, and thus complete the proof of Proposition 4.7. 

For any k large enough, let /i^ E 7Wi(M+) be the normalized restriction of /i 
to [k~^,k], so that Supp(/iA;) C M+ \ {0} is compact. The monotone convergence 
theorem yields that {^k)k converges to ;U as A; — t- 00. To approximate we have to 
stay in the class of constrained measures A^^^2(^-)i ™d thus to proceed a bit more 
carefully. To this aim, choose two sequences (aA:)^ and {bk)k satisfying Uk < < 
and 

1) Ofc decreases to inf(Supp(i/)) as /c — )• 00, 

2) bk increases to max(Supp(z^)) as k ^ 00, 

3) for any k large enough, 

i^{[ak,bk]) > {1 - k-^). (4.77) 
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Since v < a, the Radon-Nikodym theorem yields / G such that 

du{x) = f{x)dx, /(x) < — xGR_. (4.78) 

vr 

We then set the probabiUty measure 

'^"'^''^ = i^^^B^) ^^^^ " ^"')'")^['^-^^] - ^"')'^)^^' (4-79) 

whose support Supp(j^fc) C M_ \ {0} is compact. {i'k)k is easily seen to converge to 
ly as k ^ oo using monotone convergence. Moreover, it follows from (4.77)-(4.78) 
and the definition (4.79) that 

z^fc < (1 - k'^)(T. 

The fact that T^/^a: and T^i^k have finite logarithmic energy for k large enough, and 
thus that (nki^^k) satisfies Assumption 4.8, will be a consequence of (4.81)-(4.82), 
see below. 

We now prove that the sequence (/i^, i'k)k satisfies (4.76). Recall that ^{fj,, v) = 
J{T^fj,,T^i') where J is as in (4.27), namely 

JilJ-k,i^k) = II log] — - — rdT^fik{z)dT^,ij.k{w) (4.80) 
././ \z-w\ 



+ JJ [2V{z) + log\z - ^\ - log\^\)dniJk{z)dT,Uk{0 

+ jj iog^^dr,i.fc(e)dr,z.fc(c) + j \og\i\dT,Vk{i)- 

First, since S is compact, we obtain by monotone convergence 
lim [flog- — - — TdT^^ik{z)dT^iik{w) 

k->-oo J J \Z — Wl 

pk pk 2 

= //^°^|T(x)-T(y)|''^(^)^^(^) 

log , ^ . dnn{z)dnfi{w). (4.81) 

\z — w\ 

Similarly, but using moreover the metric relation (3.9), the change of variables u i— )• 
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u/{l — k ^)'^ and the inequality \u — v 

lim i I log— ^—dT^UkiCjdT^UkiC) 



hm //log ^ ■ dvk[u)dvk{v) 



k^^.lak Jak \U-V 



1 ^^TT^ViT^ 

log j j dviujdviv) 

\u — v\ 

^og—^,dT,v{i)dT,u{C). (4.82) 
I? - CI 

The same arguments moreover combined with the inequality |x — ti| > |u| for (x, u) G 
M+ X M_ yield 



^lim J I (2V(z) + log \z-C\- log \C\ynfik{z)dT,MO 





lira jj |2(x — 2-v/ax — log(l + x^)) + log |x — it| — log |d/ifc(x)(iz^fc(if) 

lim / / <^2(x-2Vax-log(l + x^)) 

+ log I (1 — k~^)^x — u\ — log \u\^diJL{x)du{u) 
|^2(x — 2\fax — log(l + x^)) + log \x — u\— log \u^dyi,{x)dv{u) 

(2V(z) + log|z-^| -log|e|)fZr,^(z)dr,z.(e) (4.83) 
After that, the continuity of T,, on <?(§-) and Lemma 4.4 provide 

hm [ logl^ldnukiO = [ log\^\dTMO- (4.84) 

Finally, (4.76) follows from (4.80)-(4.84), which completes the proof of Proposition 
4.7. □ 

4.3 Proof of Theorem 3.4 (c), (d) 

We are now in position the prove Theorem 3.4 (c), (d). The following proof follows 
closely [31, Section 2.3]. 

Proof of Theorem 3.4 (c), (d). It is enough to show that for any closed set T C 

Mi{R+) X 



limsup42log|^ivPiv((/u^,i^^) < - inf J(/x,i/), (4.85) 
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and for any open set O C A1i(M+) : 

liminf-^log|z,vIPjvf(M'^,z^'^) G o)| > - mf J{fi,u). (4.86) 

Indeed, by taking F = = x £:(M_) in (4.85) and (4.86), one obtains 

lim -^logZN = - inf J{n,u) = -J{iJ.* 

the latter quantity being finite. 

Since (4.86) has been estabhshed in Proposition 4.7, we just have to show (4.85). 
We note for convenience T^^B = {{T^fi,T^:u) : G when B C M.i{M.j^) x 

f(M_). For any closed set J" C Xi(M+) x £:(M_) we have 

Pat ((/x^, i/^) G ^) < Pjv ((T*/i^, r,z.^) G clo(r,^)) , (4.87) 

where clo(T*J^) stands for the closure of T^,F in A^i(S+) x <£"(§_). Then, since 
7Wi(S+) X <?(§_) is compact so is clo(r*J^) and, by extracting a finite covering of 
clo(r^,J^) from an appropriate covering by balls, a classical argument from LDPs 
theory (see for example the proof of [19, Theorem 4.1.11]) yields from Proposition 
4.1 that 

limsup-^logfZjvPivf(r*;U^,r*i/^) G clo(r*^)U < - inf J(/x,z/). 

(4.88) 

If (//,z^) G clo(r*J^) is such that /x({(0, 1)}) = 0, then (/i,z^) G T^F. Indeed, let 
{{T^riN ,T^^n)) be a sequence in T^F with limit (/x, i^) satisfying /i({(0,l)}) = 0. 
Since is an homeomorphism from A^i(R+) (resp. i?(M_)) to j^u G A^i(S-i-) : 
;ri({(o, 1)}) = 0} (resp. £:(§-)), this provides (r/. A) G A^i(M+) x £{M.^) such that 
(//, z^) = (T^r],T^X) and moreover the convergence of ((^/Af, Aiv))^ towards (??, A). 
Since F is closed necessarily (/i, i^) G T* J^. 

As a consequence, because J{iJ,, v) = +cxd as soon as /^({(0, 1)}) > 0, we obtain 
from the relation (4.2) 

inf J(/., u) = inf J(^, z.) = inf J(/x, i.). (4.89) 

Finally, (4.85) follows from (4.87)-(4.89). The proof of Theorem 3.4 is therefore 
complete. 

□ 
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